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Abstract 



We present the result for the one-loop squared virtual QCD corrections to the W boson pair pro- 
duction in the quark-anti-quark-annihilation channel in the limit where all kinematical invariants 
are large compared to the mass of the W boson. The infrared pole structure is in agreement with 
the prediction of Catani's general formalism for the singularities of QCD amplitudes. 



1 Introduction 



The Large Hadron Collider (LHC) is expected to have a huge impact on particle physics phe- 
nomenology. Most of the processes which will be studied at the LHC need to be calculated at least 
to next-to-leading order (NLO) in QCD whereas there are some for which a theoretical prediction 
is needed to next-next-to-leading order (NNLO). Electroweak gauge boson pair production falls 
into the latter category. One of the reasons is that the increase of the centre-of-mass energy at 
the LHC with respect to the Tevatron from 1.96 TeV to 14 TeV will result in a huge boost of the 
available data. 

The importance of hadronic W-pair production is two-fold. Firstly, it is a process which al- 
lows the measurement of the vector boson trilinear couplings and therefore a comparison with the 
Standard Model (SM) predictions. Most attempts to model New Physics, such as Supersymmetry 
and Extra-dimensions in all variations, should be able to explain any deviations by consistently ad- 
justing the anomalous couplings and/or by incorporating decays of new particles into vector boson 
pairs [1]. 

Secondly, hadronic W pair production is important for investigations of the nature of the Elec- 
troweak symmetry mechanism by contributing the dominant background for the Higgs boson me- 
diated process (see Refs. [2-14]), 

pp^H ^W*W* ^IvFV , 

in the Higgs mass range between 140 GeV < Mm < 180 GeV [15]. 

The interest in hadronic W pair production is well displayed by the fact that the Born cross 
section was calculated some thirty years ago [16]. The NLO QCD corrections were computed in 
the 90's and seen to contribute a 30% [17-21]. Next, soft gluon resummation effects were consid- 
ered in Ref. [22] whereas massless fermion-boson scattering was studied at NNLO in Ref. [23]. 
The first step towards a complete NNLO study is the computation of the NNLO two-loop virtual 
corrections in a high energy expansion, My^ ^ s, t, u [24]. The purpose of the present paper is to 
complete this study by providing the one-loop squared contributions in the same limit. 

The method used in [24], and already sketched in [26-30], is somehow different to the one used 
for the present work. In fact this time we have used the helicity matrix formalism to reduce the 
problem to a small set of integrals. We treat these again with Mellin-Bames representations [31,32] 
which are constructed by means of the MBrepresentation package [33] and then analytically 
continued in the number of space-time dimensions D = 4 — 2£ using the MB package [34] . After the 
asymptotic expansion in the mass parameter, contours are closed and integrals finally resummed 
either with the help of XSummer [35] or the PSLQ algorithm [36]. 

In Section [2l we introduce our notation, in Section [3] we verify the correctness of the infrared 
pole structure by comparing with the Catani prediction [37]. We present our results in Section |4] 
after which we conclude in Section [51 
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2 Notation 



Although the notation adopted here is identical to that of Ref. [24], we shall recapitulate it for 
completeness. The charged vector-boson production in the leading partonic scattering process 
corresponds to 

qj{Pl)+qjiP2) W-{p3,m)+W+{p4,m), (1) 

where p, denote the quark and W momenta, m is the mass of the W boson and j is a flavour index. 
We are considering down-type quark scattering in our paper. Obtaining the corresponding result 
for up-type quark scattering is actually trivial as we will show in the following. Energy-momentum 
conservation implies 

p1 + P2 = P3+P4- (2) 
We consider the scattering amplitude 9\{ for the process ([T]) at fixed values of the external parton 
momenta p,, thus p\ = P2 = ^ and = P4 = m^- The amplitude fAf may be written as a series 
expansion in the strong coupling tts. 



|fAf) 



(3) 



and we define the expansion parameter in powers of CXs(/v^) / {2tz) with /j being the renormalisation 
scale. We work in conventional dimensional regularisation, J = 4 — 2e, in the MS-scheme for the 
coupling constant renormalisation. 

We explicitly relate the bare (unrenormalised) coupling to the renormalised coupling ttg by 
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where we set the factor Se = (47r)^exp(— eye) = 1 for simplicity and P is the QCD (3-function 
known at present up to the four-loop level [38, 39] 
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The color factors in a non-Abelian SU(A^) -gauge theory are = N,Cp = {N^ — l)/2N and Tp — 
1/2. Throughout this paper, N denotes the number of colors and n^: the total number of flavors. 

For convenience, we define the function (e, m, sj^^) for the squared amplitudes summed over 
spins and colors as 



Y^\!M{qj+qj^W+ + W-)\^ = R{e,m,s,t,^i) . 
is a function of the Mandelstam variables s, t and u given by 



S={P\ +P2] 



f={Pl-P3] 



m 



{Pl -pa] 



m 



and has a perturbative expansion similar to Eq. ([3]), 



Jl(e,m,s,t,ij) 



(6) 



(7) 



(8) 
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In terms of the amplitudes the expansion coefficients in Eq. ([8]) may be expressed as 

^(0) = (i^(0)|5tf(0)), 

^(1) = ((5Vf(«)|5lf(l)) + (5lf(l)|fAf(0))) , 



(9) 
(10) 

(11) 



Note that in the above equations, only the two-loop amplitude needs to be renormalised. 
^ (0) is given by 



(0) 



A^<!ci 

+C2 
+C3 



16(l-e)^^ + 4(3-4e)± + l^ 
(1 — x) iris 

-24+ 16x+ 16e(2-x) +4^ ^ H ^ ^ ^ 



-24(1 -jc(l -x)) + 16£(2-jc(l -x)) H ^ ^ ^ ^ 



mt 



(12) 



where we have defined x = nis = ^ and only the leading physical powers (i.e. down to the 
constant) in the /M^-expansion are retained. Notice that, once the actual values of the c,- are substi- 
tuted, the terms singular in cancel as required by unitarity. This will be the case for the one-loop 
squared expression as well. The coefficients ci, C2 and C3 are in their essence combinations of EW 
coupling constants defined as 
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(13) 



The expressions for have been presented e.g. in Refs. [17, 18], whereas the real part of the 
full two-loop contributions, namely the result for the last two terms in Eq. (fTTI) in the high energy 
limit, was presented in Ref. [24] (the leading color coefficient of (f5Vf ^^^) was discussed in 
Ref. [25]). Here we provide the result for the one-loopCSone-loop result in the high energy limit, 
namely the NNLO contribution {M ^^"i) m 

In order to compute the {M (^)) we will use the helicity matrix formalism, namely we 
will express the result in terms of helicity amplitudes, M^{Xi,X2:S,t). The quark and anti-quark 
have opposite helicities in the centre-of-mass system so one helicity label above, g = ±1, suffices. 
X\ and X2 stand for the helicities of the and respectively. 
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Starting from the one-loop amplitude, the initial expression can be rearranged as 

\M^^^) = Y,Q{s,t,u)lj{s,t,u;ii,^)Mj{{pk},g), (14) 

where the C,- are coefficients, the are one-loop dimensionally regularized scalar integrals, are 
helicity matrix elements, g = ± and k= 1, ...,4. The ten helicity matrix elements Mj{pk,g) = Mj 
have been taken as defined in Ref. [40] (see also [41]): 







l), 


Mf 






Mi 


= v{p2)^i(PgU{pi)e2- P3, 




Mi 


= -viP2)^'Pguipi)ei-p4, 




Mi 


= v{p2)£/l'PgU{pi)£2-PU 




Mi 


= -v{p2)^TgU{pi)Zi-p2, 




Mi 


= v{p2)Pl3'PgU{pi)Z\-p2^2- 


P\. 


Mi 


= v{p2)pf3'PgU{p\)t\-p2^2- 


P3, 


Mi 


= v{p2)l^3'PgU{pi)Z\-pA^2- 


Pi, 


Mi 


= v{p2)Pl3'PgU{pi)Z\-pA^2- 


P3, 



where (Pg = T± = AH colour indices as well as the arguments of the polarization vectors, 

£i(p3,?ii) and e2(/'4,?i2), have been suppressed. 

Even though the representations in Eq. (15) have been used internally, we present our result 
only for the amplitude squared and summed over helicities. 



3 Infrared Pole Structure 

In the case of one-loop QCD amplitudes, their poles in £ can be expressed as a universal combi- 
nation of the tree amplitude and a colour-charge operator l(')(£). The generic form of I*^^''(£) was 
found by Catani and Seymour [42] and it was derived for the general one-loop QCD amplitude 
by integrating the real radiation graphs of the same order in perturbation series in the one-particle 
unresolved limit. 

The pole structure of our one-loop expression is given, according to the prediction by Catani, 
by acting with the operator I^^^ (£) onto the tree-level result: 

= l(i)(£)|fAf(0)) + |fAfil), (16) 

where I^^)(£) is defined as 
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The expression for the one-loop squared resuh then takes the form: 



^NNLO(lxl) ^ (18) 

= |l(i)(e)p(^ («)|^ (°)) +2Re \l(^\en^ + {^^iJj^^l) ■ 



We checked that our result is in agreement with the Catani prediction. 



4 Results 

In this section, we give explicit expressions for the finite remainder of the one-loop squared con- 
tribution 9^jjilg^^ defined as 

^^^l:^\s,t,u,m,^^) = Jl''''^^^'-'\^ (19) 

where cj^l^i^]\s,t,u,m,ij) is given by the first two terms of Eq. (fTSi) and is expanded through to 
(1), which means that it contains finite contributions as well. 

The EW structure of the finite remainder for a down-type quark can be factorised as 



^„iL^L. = A^c/ £ c,^,;i2K,^.3)- (20) 

(=1,3 

This decomposition allows one to easily obtain the result for the up-type quark scattering. The 
latter is then given by 

^i;^ = NC/l^c^jllfim,,.,^^), (21) 

1=1,3 

where one needs to use the following formulae (y = ~f ) 

n,up K'^' -2) = ri^.JM-'^^y^ -2) ' (22) 

fj. fj 

^(1x1)/ ^ \ .,(1x1) / ^ \ 

(1x1)/ ^ \ (1x1) / S 

H,up \fn,,x,-^) = HldoLi^n^.y,-^), (24) 

and naturally to make the corresponding changes in the definitions of the couplings c\, 02 and C3 
namely to use the up-type quark charge and isospin. In the following we will suppress all indices 
that indicate the type of scattered quark. 

Our result reads: 



fixn 64(l-jc)jc 192 



,211 1X4/32 3 , 
4 + l4 + 8 - + ^-1-- L 



X/ X X 1 X / \ X X 1 X 



6 _ 5 \ 2 1 1 , 1 



X 



+ 4 - + 15 - -16 ^-^ + 1- hrM L 



l—x/ \ X X^ X^ 1 —X J j ^ 



2 \ T 2 



-4 (^-32x - + 83^ + 128 (^jc - + 2^ L,„ 
^^(1,1) ^ 64( 1 - x)x ^ J_ 2 _ 2x + 3) ] + [-32(x - 2)L; - 32 f X - + 2 

-32 ( -9x--^ + 14 1 +64 (x--^ + 2 )L„ 



l—x / \ l—x 

where Lm and Ly are defined as 

Lm = Log {ms) , Ly = Log ( 1 - x) . (26) 



5 Conclusions 

We have computed the one-loop squared (a^) corrections to the process qq W+H^^ in a high 

energy expansion through to the zeroth-order in We checked that the infrared structure of our 
result agrees with the prediction of Catani's formalism for the infrared structure of QCD ampli- 
tudes. 

The present result, given as the finite remainder of the NNLO one-loop squared virtual correc- 
tions after subtraction of the structure predicted by Catani's formalism, in combination with the 
result in Ref [24] for the two-loop amplitude, completes the calculation of the virtual corrections 
to the process in the high energy limit. In a forthcoming publication, we will derive a series ex- 
pansion in the mass and integrate both results numerically similarly to what has been done for top 
quark pair production [43]. 

To complete the NNLO project one still needs to consider 2 — 3 real-virtual contributions 
and 2 — »• 4 real ones. The real-virtual corrections are known from the NLO studies on + jet 
production in Refs. [44,45]. 
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Appendix: {M^'^^ ^inite) order 8^ 

Here we present the expression for the one-loop result, {M l^^gnite) ^'^ order for down- 
type quarks. This result completes the list of the elements needed in Eq. ([8]) in order to have the 
perturbative expansion of the amplitude up to order in the high energy limit. 



i=l,3 
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^ .[-16(l-x)x] + — [-48] + 
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l-x ' V 1 



+8( -2x---^+9 I -16 fx- -^+2 )L„ 
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1-x 
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16 



1 -X 
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+8 X 
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L^ + 16 1 
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-16 I X- 

mi 



x-T^ + 4 1 +16( 1-T^ )Li2(x) 
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LjLv + 161 —X — h t I -I- i D 1 i — ) 1^12 

\ l-x/ V l-x / V l-x 



+2 U + 8 1 , 

l-x / \ l-x 

2 28 

— (l-x)x7i'^ + y(l -x)x7i;2-8(l -x)xL^-64(l -x)x + 32(l -x)xU 



+^3(12(1 -x)x- 8(1 -x)xL,)] + 
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+ 
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1 



7t---l-14x 
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69 
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l-x/ 3 V l-x 



+8 X 



l-x 
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+/71 



4 X 



1 -X 

2 



1 -X 

+ 2 ) +8(x-2)L. 
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where and Ly are defined in Eq. (|26|) and 
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s 



Logf^V (31) 



References 

[1] CDF Collaboration, Phys. Rev. Lett. 94 (2005) 21 1801; DO Collaboration, Phys. Rev. Lett. 94 (2005) 151801 

[2] M. Spira, A. Djouadi, D. Graudenz and P. M. Zerwas, Nucl. Phys. B 453 (1995) 17 

[3] S. Dawson, Nucl. Phys. B 359 (1991) 283. 

[4] R. V. Harlander and W. B. Kilgore, Phys. Rev. Lett. 88 (2002) 201801 

[5] C. Anastasiou and K. Melnikov, Nucl. Phys. B 646 (2002) 220 

[6] V. Ravindran, J. Smith and W. L. van Neerven, Nucl. Phys. B 665 (2003) 325 

[7] S. Catani, D. de Florian and M. Grazzini, JHEP 0201 (2002) 015 

[8] G. Davatz, G. Dissertori, M. Dittmar, M. Grazzini and F. Pauss, JHEP 0405 (2004) 009 

[9] C. Anastasiou, K. Melnikov and F Petriello, Phys. Rev. Lett. 93 (2004) 262002 

[10] C. Anastasiou, G. Dissertori and F Stockli, larXiv: 0707 .23731 [hep-phi . 

[11] M. Grazzini, arXiv; 0801.3232 [hep-ph] 



10 



A. Bredenstein, A. Denner, S. Dittmaier and M. M. Weber, Phys. Rev. D 74 (2006) 013004 

T. Binoth, M. Ciccolini, N. Kauer and M. Kramer, JHEP 0503 (2005) 065 | |arXiv:hep-ph/0503094) . 

T. Binoth, M. Ciccolini, N. Kauer and M. Kramer, JHEP 0612 (2006) 046 | |arXiv:hep-ph/061117(^ . 



M. Dittmar and H. K. Dreiner, Phys. Rev. D 55 (1997) 167 | |arXiv:hep-ph/9608317| 



R. W. Brown and K. O. Mikaelian, Phys. Rev. D 19 (1979) 922. 
J. Ohnemus, Phys. Rev. D 44 (1991) 1403. 
S. Frixione, Nucl. Phys. B 410 (1993) 280. 

L. J. Dixon, Z. Kunszt and A. Signer, Nucl. Phys. B 531 (1998) 3 |arXiv:hep-ph/9803250| . 



L. J. Dixon, Z. Kunszt and A. Signer, Phys. Rev. D 60 (1999) 1 14037 | |arXiv:hep-ph/9907305| . 
J. M. Campbell and R. K. ElUs, Phys. Rev. D 60 (1999) 1 13006 |arXiv:hep-ph/9905386) . 
M. Grazzini, JHEP 0601 (2006) 095 

C. Anastasiou, E. W. N. Glover and M. E. Tejeda-Yeomans, Nucl. Phys. B 629 (2002) 255 

G. Chachamis, M. Czakon and D. Eiras, arXiv:0802.4028. Thep-phl 

G. Chachamis, Acta Phys. Polon. B 38 (2007) 3563 OarXiv:0710.3035i rhep-phll. 

M. Czakon, A. Mitov and S. Moch, Phys. Lett. B651 (2007) 147, larXiv:0705.1975 [hep-ph] 

M. Czakon, A. Mitov and S. Moch, larXiv:0707.4T39l [hep-phi. 

M. Czakon, J. Gluza and T. Riemann, Phys. Rev. D71 (2005) 073009, |heFph/04 12 164| 
M. Czakon, J. Gluza and T. Riemann, Nucl. Phys. B751 (2006) l, |hep-ph/0604T0T] 



S. Actis, M. Czakon, J. Gluza and T. Riemann, Nucl. Phys. B 786 (2007) 26 
V.A. Smirnov, Phys. Lett. B460 (1999) 397, .hep-ph/9905323, 
J.B. Tausk, Phys. Lett. B469 (1999) 225, |hep-ph/9909506| 



G. Chachamis and M. Czakon, MBrepresentation .m, Unpublished 
M. Czakon, Comput. Phys. Commun. 175 (2006) 559,'hep-ph/ 051 1200| 
S. Moch and P Uwer, Comput. Phys. Commun. 174 (2006) 759, |math-ph/0508008 



H.R.P. Ferguson and D.H. Bailey, (1992), (see e.g. http://mathworld.wolfram.com/PSLQAlgorithm.htm]) l. 



S. Catani, Phys. Lett. B 427 (1998) 161 ||arXiv:hep-ph/9802439|. 



T. van Ritbergen, J.A.M. Vermaseren and S.A. Larin, Phys. Lett. B 400 (1997) 379 
M. Czakon, Nucl. Phys. B710 (2005) 485, |hep-ph/04 11261 



K. P O. Diener, B. A. Kniehl and A. Pilaftsis, Phys. Rev. D 57 (1998) 2771 
A. Denner and T. Sack, Nucl. Phys. B 306, 221 (1988). 

S. Catani and M. H. Seymour, Nucl. Phys. B 485 (1997) 291 [Erratum-ibid. B 510 (1998) 503] 
M. Czakon, larXiv:0803 . 14001 Fhep-phl 

J. M. Campbell, R. K. Ellis and G. Zanderighi, JHEP 0712 (2007) 056 
S. Dittmaier, S. KaUweit and R Uwer, Phys. Rev. Lett. 100 (2008) 062003 



11 



